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Abstract 

In liquid helium 4, many features associated to Bose statistics have been masked by the strongly 
interacting nature of the liquid. As an example of these features, we examine the shear viscosity 
of liquid helium 4 above the A point. Applying the linear-response theory to Poiseuille's formula 
for the capillary flow, the reciprocal of the shear viscosity coefficient ry is considered as a transport 
coefficient. Using the Kramers-Kronig relation, we relate a superfluid flow in a capillary with that 
in a rotating bucket, and express l/?7 in terms of the susceptibility of the system. A formula for 
the kinematic shear viscosity v is obtained which describes the influence of Bose statistics. Using 
this formula, we study the gradual fall of ^{T) from 3.7X to Tx in liquid helium 4 at 1 atm. 

PACS numbers: 67.40.-w, 67.20.+k, 67.40.Hf, 66.20.+d 
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I. INTRODUCTION 



The shear viscosity of hquid hehum 4 has been subjected to considerable experimental 
and theoretical studies. Well below the lambda temperature Tx, various excitations of a 
liquid are strictly suppressed except for phonons and rotons. Hence, it is natural to assume 
that they are responsible for the shear viscosity appearing in the damping of an oscillating 
disc, or the drag of a rotating cylinder at T < Tx. By regarding these phonons and rotons 
as a weakly interacting dilute Bose gas, the kinetic theory of gases properly describes the 
shear viscosity of liquid helium 4 at T ^ Tx ll . 

For the shear viscosity above Tx, however, the dilute-gas picture must not be applied, 
because the macroscopic condensate, which is the basis for the above picture, has not yet 
developed. Rather, we must deal with the influence of Bose statistics on the dissipation 
mechanism of a liquid. In an ordinary liquid flowing along x-direction, the shear viscosity 
causes the shear stress F^y between two adjacent layers at different velocities 

F.y = ri—, (1) 

where t] is the coefficient of shear viscosity. In the linear response theory, r/ of a sta- 
tionary flow is given by the following two-time correlation function of a tensor Jxy{t) = 



V = dt < J.y{0)J,y{t) > . (2) 

In principle, r/ of a liquid, and therefore the effect of Bose statistics on rj is obtained by 
calculating an inflnite series of the perturbation expansion of Eq.(2) with respect to the 
particle interaction f/ 3- At first sight, however, it seems to be a hopeless attempt, because 
;he dissipation in a liquid is a complicated phenomenon allowing no simple explanation p|] 
^. In liquid helium 4, many features associated to Bose statistics have been masked by the 
strongly interacting nature of the liquid. 

In this paper, we will go back to the original phenomenon which has been known from 
1938. Superfluidity was flrst discovered in a flow through a channel and a flow through a 
capillary |5|] • In an ordinary flow through a capillary (a radius a and a length L) , the velocity 
distribution under the pressure difference AP has a form such as 
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FIG. 1: A flow through a capihary. 

where r is a radius in the cylindrical coordinates (Poiseuille flow. Fig.l). In the superfluid 
phase of liquid helium 4, even when the pressure difference vanishes (AP = 0), one observes 
a non- vanishing flow {vzir^ 7^ 0), hence 77 = in Eq.(3). Instead of Eq.(2), we will regard 
Poiseuille's formula Eq.(3) as a linear-response relation. Without loss of generality, we may 
focus on a flow velocity at a single point, for example, vir = 0) on the axis of rotational 
symmetry (z-axis). We deflne a mass flow density j = pv{r = 0), and rewrite Eq.(3) as 

. = -aa ^, a = -, (4) 

where a is the conductivity of a liquid in a capillary (p is a density, and P = Pe^) In 
general, when fluid particles strongly interact to each other, a fluid swiftly responds to the 
shear stress, thereby having a small coefficient of shear viscosity 77. In Eq.(4), r] appears in 
the denominator of a. If we apply the linear-response theory to the reciprocal l/rj, we make 
a perturbation expansion of l/?7 with respect to U. In this case, an increase of U generally 
leads to an increase of l/rj and thereby a decrease ofrj . The influence of U on rj is naturally 
built into the formalism [7|. 

The change of the starting point from Eq.(2) to Eq.(4) opens a new possibility. Using the 
Kramers-Kronig relation, we will relate the frictionless capillary flow with the nonclassical 
low in a rotating bucket. This method originated in studies of electron superconductivity 
8j. Between the above two forms of superfluidity in liquid helium 4, we flnd a parallel 
relationship to that between the electrical conduction and the Meissner effect in supercon- 
ductivity. This analogy permits to express l/r/ in terms of the susceptibility of the system, 
and enables to include the effect of Bose statistics into the perturbation expansion. 

As an application of this method, we examine the gradual fall of ri{T) with decreasing 
temperature of liquid helium 4 from 3.77^^ to Tx in 1 atm. It is often said that this r){T) 



resembles r](T) of a gas than a liquid in its magnitude and in its temperature dependence. 
But this expression is somewhat misleading, because the picture of a gas has no ground in 
this regime. Rather, we must consider this gradual fall of ri(T) as a property of the liquid 
under the strong influence of Bose statistics. 

This paper is organized as follows. Section 2. A considers the shear viscosity of liquid 
helium 4 using Kramers-Kronig relation, and gives a formula for the kinematic shear viscosity 
z/(T). Section. 2. B discusses the physical reason for the suppression of the shear viscosity by 
Bose statistics. Sec. 3 examines the shear viscosity of liquid helium 4 above Tx, and Sec. 4 
makes a comparison to experiments. Section. 5 discusses some related problems. 



II. SHEAR VISCOSITY OF A BOSE LIQUID 



A. Formalism 



For liquid helium 4, we use the following hamiltonian with the repulsive interaction U 

p p,p' <? 

where $p denotes an annihilation operator of a spinless boson. 

Let us generalize Eq.(4) to the case of the oscillatory pressure as follows 

, , ^AP(uj) 

j{uj) = -a{uj)a'^^. (6) 



The conductivity spectrum a{uj) in Eq.(6) must satisfy the following sum rule 
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a{uj)du = f{a), (7) 



where /(a) is a function determined by fluid dynamics. Equation (7) is a form of the 
oscillator-strength sum rule in terms of fluid conductivity through a capillary. The Stokes 
equation gives an expression of cr{uj) and /(a) (see Appendix. A). 

At T < Ta, a frictionless flow appears in Eq.(6). Characteristic to the frictionless flow 
is the fact that in addition to the normal fluid part o"„(ti;), the conductivity spectrum criuj) 
has a sharp peak aX u = (with a very small but finite half width Ug). Hence, one obtains 

jico) = - K(c.) + AS{co)] a'^^, (8) 
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FIG. 2: The change of the conductivity spectrum (t{uj) from (a) at T > Tx to (b) at T < T^. At 
T > Tx, (t{uj) is given by the real part of Eq.(A5). At T < Tx, the sharp peak with a half width 
iOs and an area A appears around uj = 0. is a frequency at the broad peak of cr„(a;). 

where S{uj) is a simplified expression of the sharp peak at = 0, and A is an area of this 
peak Q|. Figure 2 schematically illustrates such a change of cr(to') when the system passes 
Tx- When cr^u) has a form of o"„(ci;) + yl5(ci;), r/ in Eq.(4) is given by 
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(9) 



In view of Eq.(9), the sharpness of the peak at a; = {ug — 0) leads to the disappearance 
of the shear viscosity {ri(T) ~ 0) at T < Tx- 

Let us generalize Eq. (6) so that it includes not only a mass flow responding in phase with 
AP(co') but also a mass flow responding out of phase, cr^uj) in Eq.(6) is generalized to a 
complex number cti + ia2 as follows 

,APM 



j{uj) = - [cri{uj) + io-2{uj)] a 



L 



(10) 



where cr(co') in Eq.(6) is replaced by <Ji{uj). Instead of AP(t), we will use a flctitious velocity 

v{t) in Eq.(lO). Applying the pressure gradient to a liquid is equivalent to assuming, as an 

external fleld, a velocity v(t) satisfying the equation of motion 10| 

dvjt) ^ AP(t) 
^ dt L ' ^ ^ 

When replacing AP{u)/L with —ipujviuj) in Eq.(lO), the real and imaginary part of cr(u;) 

are interchanged, 

j{uj) = p [—uja2{uj) + iujai{uj)] a'^v{uj). (12) 



Equation (12) is a linear-response formula, in which j{uj) and v{uj) forms a perturbation 
energy such as j^u) ■ v{uj) [11]. The susceptibility x{^) consists of a real part —uja2{uj) for 
the non-dissipative flow and an imaginary part iujai{uj) for the dissipative one. 

Causality requires that one observes a mass flow only after the pressure is applied. We 
obtain the following Kramers-Kronig relation for —uja2{uj) and ujai{uj) in Eq.(12) 

Mu^) = -l ^-^73^' (13) 

u'a^i.') ^ r d.^, (14) 

If one determines o"2(ix;) in the non-dissipative flow, one obtains cri(ti;) using Eq.(13). 

As a non-dissipative flow, we consider a flow in a rotating bucket. In fluid dynamics, the 
viscous dissipation in an incompressible fluid is estimated by the dissipation function 
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$(r) = li] - -efcfc^ijj , (15) 

{tij = {dvi/dxj + dvj /dxi)/2 is the shear velocity). Different flow patterns have different 
degrees of dissipation: For the Poiseuille flow Eq.(3), the dissipation function Eq.(15) is 
not zero at every r except for r = 0. Fluid particles in a capillary flow experience thermal 
dissipation not only at the boundary, but also inside of the flow. On the other hand, in a 
rotating bucket, a liquid makes the rigid-body rotation owing to its viscosity. Except at the 
boundary to the wall, there is no frictional force within a liquid, and the flow is therefore a 
non-dissipative one. (The velocity is a product of the radius and the rotational velocity Q 
such as Vd{r) = fl x r, for which Eq.(15) is zero at every r.) 

The flow in a rotating bucket is formulated using the generalized susceptibility of the 



system x(r,a;) [121]. When we assume Vd{r,t) = fl x re*"^* as a v{uj) in Eq.(12), j{r,uj) and 
Vd{r, uj) form the perturbation energy like J (Ixv^- j- Since the flow in a rotating bucket is a 
non-dissipative system, j{r, uj) is a dynamical response of a liquid to the mechanical external 
fleld Vd{r, u), such as j{r, u) = x(r, co')vd(r, u). The influence of the wall motion propagates 
from the wall to the center along the radial direction, which is perpendicular to the particle 
motion driven by rotation. Hence, the flow in a rotating bucket is a transverse response 
described by the transverse susceptibihty x^{q,^)- In the right-hand side of Eq.(12), the 
real part —puja2{oj)a? is expressed as 



pa uja2{oj) = limx {q,uj)- (16) 

9^0 



Using Eq.(16) in the right-hand side of Eq.(13), we obtain ai{uj) for the capillary flow 

2 / A 2 /""^ limg_oX^(?,^) 
pa aiiu) = / du ^— . (17) 



Equation(17) has the following meaning. Quantum mechanics states that, in the decay 
from an excited state with an energy level E to a ground state with Eq, the higher excitation 
energy E causes the shorter relaxation time r. The time-dependent perturbation theory says 

-~|E-Eo|. (18) 
r 

In Eq.(17), the left-hand side includes the relaxation time r as di = p/(4?7) = p/^AGr) (see 
Appendix.B), whereas the right-hand side includes the excitation spectrum in x"^(g, cj). In 
this sense, Eq.(17) is a many-body theoretical expression of Eq.(18). 

In the normal fluid phase, x^iQ^^) = X^iQy^) is satisfied at small q and u, and one can 
replace ^-^(g, cj) in Eq.(17) by x-^(g, cj) for a small q and u) [l3j. Hence, the conductivity 
(Tin{uj) of the capillary flow in the normal fluid phase is given by 

pa ain{uj) = / duj ^ . (19) 



In the superfluid phase, under the strong influence of Bose statistics, the condition of 
X^{q,uj) = x^{q,uj) at g — > is violated (see Sec.2.B). Consequently, one cannot replace 
X^{q,u) by x^(g, c<j) in Eq.(17). In addition to ainioj), one must separately consider the 
contribution of x^{q,uj) — x'^(g, cj). Hence, one obtains 

/ /N / /A , 2 limg^o[x^(g,^) -x^(g,^)] .„„x 

oi[uj ) = ain[uj ) ^ -I dw , (20) 

piia J Q uj — uj 

in which limq^o[x^(Q') ^) " X^(Q')^)] corresponds to a superfluid component. In general, 

superfluidity is not rigid against external perturbations oscillating at not low frequencies: 

Above a certain cuq, limq^o[x^(?! i^) ~ X"^(Q'!'^)] becomes zero. One can see (J\{oj) of such 

a system using a simple example. We assume a constant limg^o[x'^(Q') 0) — ^^(^jO)] for 

< < cuo, and for ujq < oj. In this case, we obtain 

ai(u;) =(Ti„(u;) + ^lim[x^(g,0)-X^(g,0)]D,„(u;). (21) 



where 
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^-oM = 7r-lnK-l|. (22) 
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FIG. 3: D^^{u) for loq = 1. 



Figure. 3 shows D^J^{uJ) for ujq = 1. For every ujq, D^q{u!) resembles the (5(u;)-function. 
This simple example suggests that, for all probable forms of limg^olx^ {q,u!) — x^{q,^)], 
the sharp peak around a; = is a general feature of a"i(a;). On comparing Eqs.(21) and 
(22) with Eq.(8), we obtain an expected form of a{u;) for the superfluid flow, and A — 
21img_^o[x^(Q'; 0) — x^l?, 0)]/(p7ra^). The half width a;^ of the peak in ai{Lu) is determined 
by Uq in ~ X^j it gives us a measure for the rigidity of superfluidity. 

When the system passes Tx, cr(a;) changes from Fig.2(a) to 2(b) under the sum rule 
Eq.(7). The area A of the sharp peak is equal to that of the shaded region in Fig.2(b). 
At a; = 0, the conductivity is enhanced by the sharp peak as A/cUs, but the original a{0) 
becomes a{0) — Act by the sum rule. We approximate the shaded region in Fig. 2(b) by 
a triangle. (A broad peak of cr„(a;) is located at ujc)- We take into account the sum rule 
as —AauJc/2 + A = 0, hence Acr = 2A/uJc. After all, the conductivity at = becomes 
[cr(0) — 2A/uJc] + Ajijjs at T < T\. Using this form of cr(0) and the expression of ^4, one 
obtains r] — p/(4(7(0)) as 



Here, we define the mechanical superfluid density Ps(T) = \img^o[x^{q,0) — X^(?)0)]) 
which does not always agree with the conventional thermodynamical superfluid density 
ps{T). (By "thermodynamical", we imply the quantity that remains finite in the V ^ oo 
limit.) Using ujg <S uJc in Eq.(23), we obtain the following formula of the kinematic shear 
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(23) 
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lim[x^(g,0)-x^(g, 0)] 
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viscosity viT) = TjiT)/ p{T) 

where z/„ is the kinematic shear viscosity in the normal fluid phase, and it satisfies c'"i„(0) = 
p/{Arjn) = l/(4z/„) in Eq.(23). (The total density p in Ec[.(23) slightly and monotonically 
increases with decreasing temperature from 4.2K to Tx 1^. Comparing Eq.(24) with (23), 
we see that z^(T) = ri{T)/p{T) is a more appropriate quantity than ri{T) for describing the 
change of the system around Tx.) 

One notes the following features in the formula (24). 

(1) i^iT) of a superfiuid is expressed as an infinite power series of z/„, and the infiuence of 
Bose statistics appears in its coefficients. This result does not depend on a particular model 
of a liquid, but on the general argument. (The microscopic derivation of Un depends on the 
model, which is a subject of the liquid theory and beyond the scope of this paper.) 

(2) Because of Ug — 0, which characterizes the sharp peak in Fig.2(b), a small change 
of Ps{T) in the denominator of the right-hand side is strongly enhanced to an observable 
change of z^(T). 

(3) The existence of 1/a^ in front of psiT)/ p indicates that a frictionless superfiuid fiow 
appears only in a narrow capillary with a small radius a: A narrower capillary shows a 
clearer evidence of a frictionless fiow. 



B. Physical interpretation 

There is a physical reason to expect that the shear viscosity of a Bose liquid falls off at 
low temperature. For the shear viscosity of a classical liquid. Maxwell obtained a simple 
formula rj = Gt (Maxwell's relation) using a physical argument 1^. {G is the modulus of 
rigidity, r is a relaxation time of the process in which the fluid motion relaxes the shear 
stress between adjacent layers in a flow. See Appendix.B). In the vicinity of Tx, no apparent 
structural transformation is observed in liquid helium 4. Hence, G may be a constant at the 
flrst approximation, and therefore the fall of the shear viscosity is attributed to a decrease of 
T. In view of Eq.(18), the decrease of r suggests that the excitation energy E increases owing 
to Bose statistics. The relationship between the excitation energy and Bose statistics dates 
back to Feynman's argument on the scarcity of the low-energy excitation in liquid helium 



16| . in which he explained how Bose statistics affects the many-body wave function in 



configuration space. To the shear viscosity, we will apply his explanation. 

Consider a flow in Fig.l. White circles represent an initial distribution of fluid particles. 
The long thin arrows move white circles on a solid straight line to black circles on a one- 
point-dotted-line curve. (A viscous liquid shows such a spatial gradient of the velocity. 
The influence of adjacent layers in a flow propagates along a direction perpendicular to 
the particle motion. Hence, the excitation caused by the shear viscosity is a transverse 
excitation.) Let us assume that a liquid in Fig.l is in the Bose-Einstein Condensation 
(BEC) phase, and the many-body wave function has permutation symmetry everywhere 
in a capillary. At first sight, these displacements by long arrows seem to be a large-scale 
configuration change, but they are reproduced by a set of slight displacements by short 
thick arrows from the initial particles. In contrast with the longitudinal displacement, the 
transverse displacement does not change the particle density in the large scale, and therefore, 
to any given particle after displacement, it is always possible to find a particle being close 
to that particle in the initial configuration 17|]. In Bose statistics, owing to permutation 
symmetry, one cannot distinguish between two types of particles after displacement, one 
moved from the neighboring position by the short arrow, and the other moved from distant 
initial positions by the long arrow. Even if the displacement made by the long arrows is a 
large displacement in classical statistics, it is only a slight displacement by the short arrows 
in Bose statistics . 

Let us imagine this situation in 3 N- dimensional configuration space. The excited state 
made of slight displacements, which is a characteristic of Bose statistics, lies in a small 
distance from the ground state in configuration space. Since the wave function of the 
excited state is orthogonal to that of the ground state in the integral over configurations, 
the amplitude of the many-body wave function of the excited state must spatially oscillate 
around zero. Accordingly, it oscillates within a small distance in configuration space. The 
kinetic energy of the system is determined by the 3 N- dimensional gradient of the many-body 
wave function in configuration space, and therefore this steep rise and fall of the amplitude 
raises the excitation energy of this wave function. The relaxation from such an excited state 
is a rapid process with a small r. This mechanism explains why Bose statistics leads to the 
small coefficient of shear viscosity rj = Gt. 

When the system is at high temperature, the coherent wave function has a microscopic 



10 



3 



^ 2 



£ 1 







2 

T(K) 



FIG. 4: The temperature dependence of the kinematic shear viscosity J^(T) of liquid hehum 4 



obtained in Ref. 



size. If a long arrow in Fig.l takes a particle to a position beyond the coherent wave function 
including that particle, one cannot regard the particle after displacement as an equivalent of 
the initial one. The mechanism below Tx does not work for the large displacement extending 
over two different wave functions. Hence, the relaxation time r changes to an ordinary long 
r, which is characteristic of a normal liquid [l^. 



III. SHEAR VISCOSITY ABOVE Tx 



Figure 4 shows i^iT) of liquid helium 4 in the vicinity of Tx [19]. One notes that z^(T) 
does not abruptly drops to zero at Tx- Rather, after reaching a maximum value at 3.7K, it 
gradually decreases with decreasing temperature, and finally drops to zero at the A-point. 
z/(T) of a classical liquid in Eq.(24)) has the general property of increasing monotonically 
with decreasing temperature 20|]. In Fig. 4, u above 3.7 K seems to follow this pro per ty, 
whereas its gradual fall below 3.77^ suggests Ps(T) 7^ in Eq.(24) at Tx < T < 3.7K \21 \. 

In an ideal Bose gas, one knows Ps{T) = at T > Ta, hence obtains z^(T) = z/„ at 
T > Ta in Eq.(24). This means that, without the interaction between particles, BEG is the 
necessary condition for the superfluid flow. To explain the gradual fall of z^(T) just above 
Ta, we must obtain x^il,^) ~X^{1j^) under the repulsive interaction U. The susceptibility 
is decomposed into the longitudinal and transverse part (/x, u = x,y, z) 



r 



(25) 
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FIG. 5: Exchange of particles in (a) between J^Jjy (a large bubble) and a bubble excitation (an 
inner small one) yields (b). 

For the later use, we define a term proportional to g^gj, in Xfiv by Xf_^y = {q^iqu / 'i^)[x^ {q, ^^) — 
X^(g, cj)] for considering Ps(T) = \\mq_Q[q^ / {q^qy)\x^v Under the repulsive interaction 
Hi{r), Xfj_y is derived from 



<G\T^J^{x,t)J,{QMG>-- 



< Tt- J^(x, r) J,^(0, 0)exp 


i-f 

. Jo 




|0 > 


< exp 


- / dTHi{T) 
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|0 > 



(26) 



where J^(g,r) = Ep,niP + Q/'^)>^H^p+,(^~''^"" ^ (/^ = V(^b2^),^ = ^t). Compared to 
an ordinary liquid, liquid helium 4 above Tx has a 10^^ times smaller coefficient of shear 
viscosity as shown in Fig.4. Although in the normal phase, it is already an anomalous 
liquid under the strong infiuence of Bose statistics. Hence, it seems natural to assume that, 
as T —>■ Tx in the normal phase, large but not yet macroscopic coherent wave functions 
gradually grow as a thermal equilibrium state, and suppress the shear viscosity j22| . 



Here, we summarize the result of Ref. 23j, which formulates the interpretation in Sec.2B. 
Figure 5 illustrates the current-current response tensor J^(x, t)J^{0, 0) (a large bubble with 
/i and u) in the liquid: Owing to exp(— / Hi{T)dr), scattering of particles frequently occurs 
as illustrated by an inner small bubble with a dotted line U. As the order of the perturbation 
increases, J^(x, r) Jj,(0, 0) in the vicinity of Tx gradually reveals a new effect due to Bose 



statistics: The large bubb 
wave function whole 



e J/i(x, r)X(0, 0) and the small bubble in Fig. 5 (a) form a coherent 



24l |: When one of the two particles in the large bubble and in the 



small bubble have the same momentum [p = p'), and when the other in both bubbles have 
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another same momentum {p+q — p'+q'), Bose statistics forces us to include in the expansion 
a graph made by exchanging these particle. This exchange yields Fig. 5(b), in which two 
bubbles with the same momenta are linked by the repulsive interaction. With decreasing 
temperature, the coherent wave function grows, and such an exchange of particles occurs 
many times. Furthermore, among possible physical processes, processes including p = 
particles get to play a more dominant role than others: In Fig.5, a bubble with p — 
corresponds to an excitation from the rest particle, and a bubble with p — —q corresponds 
to a decay into the rest one. Taking only such processes and continuing these exchanges, 
one obtains x^i/(5, 0) for Ps{T) 

x.'^M-^y ,_UF,iqy (27) 

where 

p . ) ^ (exp(/3[S - /x]) - 1)-^ - (exp(/3[e(g) + S - /x]) - 1)"^ 

is a positive and monotonically decreasing function of q^, which approaches zero as — > oo. 
An expansion of -F/3(g) around = 0, Ff){0) — bq'^ + ■ ■ • has a form such as 



4sinh^(MMn^ 



P 1 q' , 

^anh('^tMT)-E]|^2^ 



(29) 



(a) With decreasing temperature, /i gradually approaches E. As // — E — > 0, F/3(?) in 
Eq.(28) increases, and it creates a divergence in Eq.(27): This divergence first occurs at 

q^ = 0, because Fp{q) is a positive decreasing function of q^. For a small g, 1 — UFplq) in 
Eq.(27) is approximated as [1 — ^/^^(O)] + Ubq^. As /x — E — >• 0, UFp{Q) increases and finally 
reaches 1, that is, 

(;^ = 4rf(*m^). (30) 

At this point, the denominator in Eq.(27) gets to begin with g^, and Xn^iQ^ 0) therefore has 
a form of q^q^/q^ at g ^ 0, the coefficient of which is Fpi^Q)/ {2VUh). By the definition of 
Ps{T) = lvi^q^Q{.q^ / qnqu)Xnv, one obtains 

„ 1 m 

^^(^ " T^sinh|/3[//(r)-E]r ^^'^> 
with the aid of Eqs.(29) and (30). Here, we call T satisfying Eq.(30) the onset temperature 
Ton of the nonclassical behavior. 
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(b) In the vicinity of Tx, Eq.(30) is approximated as UP = /3^[/i(T) — for a small 
fi — Ti. This condition has two solutions yu(T) = S(f/) ± ^/UksT. The repulsive Bose system 
is generally assumed to undergo BEC as well as a free Bose gas. Hence, with decreasing 
temperature, /i(T) of repulsive Bose system should reach S(?7) at a finite temperature, 
during which course the system necessarily passes a state satisfying yu(T) = S(f/) — y/UkBT. 
Consequently, Ton is always above Tx, and i^(T) deviates from z/„ just above Tx in Eq.(24). 

(c) As T approaches Tx from Ton, Ps(T) approaches (m/\^) [exp(/9[/i(T) — S]) — l]^"*^. This 
means that at T = T^, Ps{Tx) agrees with the conventional thermodynamical superfluid 
density ps{Tx), and it abruptly grows to a macroscopic number. 



(d) At Tx < T < Ton, Eq.(31) serves as an interpolation formula for Ps(T) [25|. 



IV. COMPARISON TO EXPERIMENTS 



At an early sta ge i n the study of superfluidity, ri{T) was measured by experiments using 
the capillary flow 26 1 28]. But early measurements have been superseded by that using 



more accurate techniques such as a vibrating wire [2^ 3^ [3l|. Above the A point, these 
different experimental methods give us quantitatively similar data on 1]{T). Using these 
accumulated data, a precise curve of ri{T) and z^(T) above Tx was obtained by statistical 



analysis such as cubic spline fits 19l- We will compare our formula with the results in Ref 



[1) In Eq.(24), we assume Ton = 3.7-?^ and Un = ^(Ton)- We obtain p^iT)/ p by 



Ps{T) 



u{T)-^ - u{To, 



(32) 



In the rotation experiment by Hess and Fairbank 32], the moment of inertia Iz just above 
Tx is slightly smaller than the normal phase value If (see Sec. 5. A). Using these data, Ref. 23] 
roughly estimates ps{T)/p as p,(TA + 0.03ip/p = 8 x 10-^ and ps{Tx + 0-28K)/ p ^ 3 x 10"^ 
The previous experiments in Ref. 26| 27|] 28|] differ in the capillary radius a. If we use a 
typical value a ~ 5 x 10~^cm, with these pg{T)/p in Ref 23| and z^(T) in Fig.4, we obtain a 
rough estimate of cu^ as cu^ ^ 4 x 10~^rad/s. Figure 6 shows the temperature dependence of 
Pg(T)/p derived by using z/(T) (Fig.4) in Eq.(32). (The absolute value is adjusted to p^/p 
at Tx + 0.037^ and Tx + 0.28K in Ref |23j]. Since the precision of these values derived from 
currently available data is limited, the absolute value of Pg{T)/p in Fig.6 has a statistical 
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FIG. 6: ps{T)/p obtained by Eq.(32) using v{T) in Fig.4. 



uncertainty.) 

(2) In Sec. 3, we obtained the interpolation formula for Pg{T)/p, Eq.(31). Here, we assume 
— E changes with temperature according to the formula 



{ga{x) — '^n^^ /n"'), where the particle interaction U and the particle density p of liquid 
hehum 4 are renormalized to Tx — 2.17 K. The temperature dependence of Pg{T)/p derived 
from Eqs.(31) and (33) bears a qualitative resemblance to pg(T)/p in Fig.6. But, although 
remaining very small at T\ + 0.2K < T < 3.7K, it abruptly grows at 2.3K, and reaches a 
macroscopic number at Tx (it resembles the shape of the letter L) . Experimentally, however, 
as temperature decreases from 3.7K to Tx, Pg{T)/p gradually grows as in Fig.6. This means 
that Eq.(31) is too simple to compare with the real system. With decreasing temperature, 
in addition to the particle with p — 0, other particles having small but finite momenta get to 
contribute to the divergence of x^uiq,^) as well. (In addition to Eq.(28), a new F^lq) 
including p 7^ also satisfies 1 — t/F^(0) = in Eq.(27).) The participation of p 7^ particles 
into Ps{T) is a physically natural phenomenon. For the repulsive Bose system, particles are 
likely to spread uniformly in coordinate space due to the repulsive force. This feature makes 
the particles with p ^ behave similarly with other particles, especially with the particle 
having zero momentum. // they move at different velocities along the fl,ow direction, the 




3/2 



2 



- 1 



(33) 
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particle density becomes locally high, thus raising the interaction energy. This is a reason 



why many particles with p ^ participate in a superfluid even above T\ 33 1. 

(3) For the rigidity of superfluidity, we used two types of phenomenological parameters 
in Sec. 2. Among them, ujs of the sharp peak in cri(co') is useful for comparing the theory 
with experiments. The quantity which has a physically clear meaning, however, is ujq in 
\\mq^Q[x^{q,uj) — c<j)]. As an example, we used a step function of u with a width 
ujq. Its Hilbert transform cri(co') in Eq.(22) shows a 5-function-like peak with a width uJs 
such as UJs = lO^^c^o ~ in Fig.3. Since Us = 10~^rad/s, the critical frequency ujq 
above which the condition of superfluidity limg^o i-i-') 7^ lini^^o i^) is violated is 
lO^^rad/s ~ Irad/s. 

(4) Equations (30) and (33) with Ton = 3.7/^ gives us a rough estimate of f/ as f/ ~ 
3.4 X 10~^^erg, which is somewhat larger than the value obtained from the sound velocity 
using the Bogoliubov formula. 

V. DISCUSSION 

A. Superfluidity in the dissipative and the non-dissipative systems 

On the onset mechanism of superfluidity, one can see physical differences between dis- 
sipative and non-dissipative flows. As an example of non-dissipative flows, we considered 
the flow in a rotating bucket in Sec. 2, in which a quantity directly indicating the onset of 
superfluidity is the moment of inertia 

1 - (34) 



p 

where If is its classical value. Equation (34) and (24) have the following differences. 

(a) In Eq.(34), Ps{T) appears as a correction to the coefficient of the linear term If, 
whereas Eq.(24) has a form of infinite power series of Vn and Ps{T)- With decreasing tem- 
perature, the higher-order terms become dominant in Eq.(24). 

(b) In Eq.(34), the change of ~ affect s L without being enhanced, and therefore 



the small Ps(T) only slightly affects above Tx |23|. In Eq.(24), the change of does 
not directly affect u, but through the dispersion integral as in Eq.(20). By this mechanism, 
a small change of the system leading to superfluidity is enhanced to an easily observable 
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scale, which amphfication mechanism does not appear in superfluidity in the non-dissipative 
systems. 

(c) Equation (34) includes no phenomenological parameter, whereas Eq.(24) includes 
the parameter representing the measure of the rigidity of superfluidity. This feature is 
characteristic of superfluidity in the dissipative systems. The stability of a superfluid has 
been studied in the context of the critical velocity of the persistent current js^. To obtain 
ujs or tuo, similar considerations are needed for the stability of a superfluid against oscillating 
perturbations such as the oscillating velocity f (r)e*'^*. This is a future problem. 

Anomalously high thermal conductivity of liquid helium 4 at T < Ta is another example 
of superfluidity in the dissipative system. Heat flow q is expressed as g = — /tVT, where k 
is the coefficient of thermal conductivity. In the critical region above T\ {T /T\ — 1 < 10^'^), 
the rapid rise of k is observed. At T = Tx, k jumps to an at least 10'' times higher value 
than K just above T\. In Tx < T < 3.7 K, however, k does not show a symptom of its 
rise, in contrast with r] which shows a symptom of its fall in the same region js^]. In the 
case of the capillary flow, we can assume the flow in a rotating bucket as a non-dissipative 
counterpart, whereas in the case of heat conduction, we do not know such a counterpart. 
Theoretically, the coefficient of thermal conductivity is expressed by a correlation function 
which has a formally similar form to that of the shear viscosity. But heat conduction 
is a transport phenomenon of energy. The scalar fleld (temperature field) has only a small 
variety of spatial distribution compared to the vector field (velocity field), and therefore heat 
conduction is always a dissipative phenomenon. Hence, we can not apply the formalism in 
this paper to the onset mechanism of the anomalous thermal conductivity. This formal 
difference between the shear viscosity and the thermal conductivity is consistent with the 
experimental difference between them at T\ < T < 3.7 K. 

B. Comparison to Fermi liquids 

In liquid helium 3, the fall of the shear viscosity at is known as a parallel phenomenon 
to that of liquid helium 4. The formalism in Sec. 2 is applicable to the shear viscosity of 
liquid helium 3 as well. For the behavior above Tc, however, there is a striking difference 
between liquid helium 3 and 4. The phenomenon occurring in fermions in the vicinity of 
is not a gradual growth of the coherent wave function, but a formation of the Cooper pairs 
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from two fermions. (This difference evidently appears in tlie temperature dependence of the 
specific heat: C{T) of hquid hehum 3 above Tc does not show a symptom of its rise.) Once 
the Cooper pairs are formed, they are composite bosons situated at low temperature and 
high density, and immediately jumps to the ESP or BW state. Hence, the shear viscosity of 
liquid helium 3 shows an abrupt drop at without a gradual fall above T^. 

In electron superconductivity, the fluctuation-enhanced conductivity cr' is observed above 
Tc (see the text by M.Tinkham in [s!]). In bulk superconductors, however, the ratio of a' 
to the normal conductivity cr„ is about 10~^ at the critical region, and zero outside of this 
region. Practically, it is unlikely that thermal fluctuations create a large change of a at 
temperatures outside of the critical region. 



C. Velocimetry technique 

Since the discovery of superfluidity in liquid helium 4, a lot of measurements have been 
done, but there is so little direct experimental information about the flow patterns. Among 
various visualization techniques used in ordinary liquids, particle image velocimetry (PIV), 
which records the motion of micrometre-scale solid particles suspended in the fluid as tracer 
particles, recently becomes available in superfluid helium 4. Until now, this technique has 



been mainly used for the study of turbulent flows [36| [37| , but it has the potential to give 
us more information. The quantity which directly indicates the onset of superfluidity in 
dissipative system is the change of the conductivity spectrum a{uj). In view of the gradual 
decrease of r] above Tx, the sharp peak at u; = and the corresponding change from cr{uj) 
to Cniuj) in Fig. 2 must already appear at T > Ta. To confirm this prediction, under the 
slowly oscillating pressure, time- resolved measurement of the oscillating flow velocity v{r, t) 
is necessary. Such an experiment must be performed in a thin capillary with an inner radius 
of 10~^ ~ 10"^ mm. If the PIV experiment is performed under such conditions and the 
data of v{r = 0,t) is Fourier-transformed, the change of cr{u!) in Fig.2 will be observed. 



Since quantized vortices are absent in the capil 



the normal fluid part and trace its velocity 



ary flow, tracer particles interact only with 



381]. Hence, in Fig.2 (b) only the change from 



cr{uj) to (Jniuj) will be observed. The quantity A in Eq.(8) is equal to the area of the shaded 
region cr{u) — (j„(a;) in Fig. 2(b). Using thus obtained A, one can determine ujg by comparing 
the experimental data of f]{T) with Eq.(9). Furthermore, using ira'^A = 2p^{T)/p, one will 
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obtain a new estimate of Pg{T)/ p aX, T > T\. Such an experiment may be a difficult one, 
but it will provide us valuable information. 



APPENDIX A: CONDUCTIVITY SPECTRUM a{uj) 

The Stokes equation under the oscillatory pressure gradient APe*"^*/-^ written in the 
cylindrical polar coordinate as follows 



dv _ f d d \ APe^'^* 
dt \dr'^ rdr J pL 

The velocity has the following form 

APe*'^* 

'v{r,t) = — — h Av(r,t), 

lUjpL 

with the boundary condition of v(a, t) = 0. Av(r, t) satisfies 

dAv{r,t) _ ( d ^ \ /\ ( \ 
dt V dr'^ rdr ) ' ' 



(Al) 



(A2) 



(A3) 



and therefore Av(r, i) has a solution written by the Bessel function Jo(iAr) with A = 
{l + ^)^JuJ/{2v). Hence, 



v{r,t) = 



APe 



iuit 



1 - 



Jo{i\r] 



(A4) 



iujpL \ jQ[iXa) 

At r = 0, the conductivity spectrum cr(a;), which is given by pv{0,t) = cr(a;)a^APe*'^*/-^) 
has the following form 

/ \ 



1 - 



v 



Jo 




(A5) 



The real part of Eq.(A5) gives a curve of a{io) in Fig.2(a) (Re (j(0) = p/ (4?7)), and determines 
/(a) cx in Eq.(7). {u disappears in J a{ijj)du}.) Im cr(a;) gives an expression of (J2{uj) 
in Eq.(lO), but gives no concrete form of 0-2(^7) in Eq.(16), because it is derived from the 
phenomenological equation with dissipation like the Stokes equation. 

APPENDIX B: MAXWELL'S RELATION 

Consider shear transformation of a solid and of a liquid. In a solid, shear stress F^y is 
proportional to a shear angle as F^y — G(f), where G is the modulus of rigidity. The value 
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FIG. 7: In a liquid flowing along the rc-direction, owing to the velocity gradient along the y- 
direction, a small rectangular part of a liquid is deformed to a parallelogram. 

of G is determined by dynamical processes in which vacancies in a solid move to neighboring 
positions over the energy barriers. As (j) increases, F^y increases as follows, 

^ = G^. (Bl) 

In a liquid, the motion of a fluid changes the relative positions of particles, reducing the 
shear stress F^y to a certain value. Presumably, the rate of such relaxations is proportional 
to the magnitude of F^y, and one obtains 

dFxy ^d(f) Fxy /T30^ 

-df^^Tt-— ^^^^ 

where r is a relaxation time. In the stationary flow after relaxation, F^y remains constant, 
and one obtains 

= (B3) 
dt T ^ ' 

Figure. 7 shows two particles 1 and 2, each of which starts at {x^y) and (x, y + Ay) simul- 
taneously and moves along the a;-direction. Assume that there is velocity gradient Vx{y) 
along y direction. After At has passed, they (1' and 2') are at a distance of Au^-At along 
the x-direction. The shear angle increases from zero to A0 as a result, which satisfies 
Ai>3;Ai = AyA0 as depicted in Fig. 7. Hence, we obtain 

dy dt' ^ ' 

Substituting Eq.(B4) into Eq.(B3), and comparing it with Eq.(l), one obtains r] = Gr 
(Maxwell's relation). 
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